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Abstract. We show that the n-fold integrals x of the magnetic susceptibility 
of the Ising model, as well as various other n-fold integrals of the "Ising class" , 
or n-fold integrals from enumerative combinatorics, like lattice Green functions, 
correspond to a distinguished class of function generalising algebraic functions: 
they are actually diagonals of rational functions. As a consequence, the power 
series expansions of the, analytic at x = 0, solutions of these linear differential 
equations "Derived From Geometry" are globally bounded, which means that, 
after just one rescaling of the expansion variable, they can be cast into series 
expansions with integer coefficients. We also give several results showing that the 
unique analytical solution of Calabi-Yau ODEs, and more generally of MUM 
ODEs, are almost always diagonal of rational functions. Besides, in a more 
enumerative combinatorics context, generating functions whose coefficients are 
expressed in terms of nested sums of products of binomial terms can also be 
shown to be diagonals of rational functions. We finally address the question of 
the relations between the notion of integrality (series with integer coefficients, or, 
more generally, globally bounded series) and the modularity of ODEs. 
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1. Introduction 

The series expansions of many magnetic susceptibilities (or many other quantities, 
like the spontaneous magnetisation) of the Ising model on various lattices in arbitrary 
dimensions are actually series with integer coefficients [1, 2, 3]. This is a consequence 
of the fact that, in a van der Waerden type expansion of the susceptibility, all 
the contributing graphs are the ones with exactly two odd-degree vertices and the 
number of such graphs is an integer. When series expansions in theoretical physics, 
or mathematical physics, do not have such an obvious counting interpretation, the 
puzzling emergence of series with integer coefficients is a strong indication that some 
fundamental structure, symmetry, concept have been overlooked, and that a deeper 
understanding of the problem remains to be discovered^. Algebraic functions are 
known to produce series with integer coefficients. Eisenstein's theorem [5] states 
that the Taylor series of a (branch of an) algebraic function can be recast into a 
series with integer coefficients, up to a rescaling by a constant (Eiscnstcin constant). 
An intriguing result due to Fatou [6] (see pp. 368-373) states that a power series 
with integer coefficients and radius of convergence (at least) one, is either rational, 
or transcendental. This result also appears in Polya and Szego's famous Aufgaben 
book [7] (see Problem VIII-167). Polya [8] conjectured a stronger result, namely 
that a power series with integer coefficients which converges in the open unit disk is 
either rational, or admits the unit circle as a natural boundary (i.e. it has no analytic 
continuation beyond the unit disk). This was eventually proved^] by Carlson [10]. 
Along this natural boundary line, it is worth recalling [11, 12, 13, 14, 15] that the 
series expansions of the full magnetic susceptibility of the 2D Ising model corresponds 
to a power series with integer coefficients]". For them, the unit circle certainly arises as 
a natural boundary [16] (with respect to the modulus variable k), but, unfortunately, 
this cannot be justified by Carlson's theoremft. 

A series with natural boundaries cannot be D-finite^ i.e. solution of a linear 
differential equation with polynomial coefficients [20, 21]. For simplicity, let us restrict 
to series with integer coefficients (or series that have integer coefficients up to a 
variable rescaling), that are series expansions of D- finite functions. Wu, McCoy, 
Tracy and Barouch [22] have shown that the previous full magnetic susceptibility 
of the 2D Ising model can be expressed (up to a normalisation factor (1 — s) 1 / 4 /s, 
see [13, 23]) as an infinite sum of n-fold integrals, denoted by x i which are actually 
D-finite^,. We found out that the corresponding (minimal order) differential operators 
are Fuchsian [11, 13], and, in fact, "special" Fuchsian operators: the critical exponents 

| The emergence of positive integer coefficients corresponds to the existence of some underlying 
measure [4]. 

% The Polya-Carlson result can be used to prove that some integer sequences, such as the sequence 
of prime numbers (pn) [9], do not satisfy any linear recurrence relation with polynomial coefficients, 
f In some variable w [12, 11, 13, 14]. In the modulus variable k, one needs to perform a simple 
rescaling by a factor 2 or 4 according to the type of (high, or low temperature) expansions, 
ft The radius of convergence is 1 with respect to the modulus variable k, in which the series does not 
have integer coefficients, being globally bounded only (this means that it can be recast into a series 
with integer coefficients by one rescaling of the variable k). If one considers the series expansion with 
respect to another variable (such as w) in which the series does have integer coefficients, then the 
radius of convergence is not 1. 

ft D-finite series are sometimes called holonomic. A priori, these notions differ [17]. The equivalence 
of these notions is proved by profound results of Bernstein [18] and Kashiwara [19]. 
§ For Ising models on higher dimensional lattices [1, 2, 3] no such decomposition of susceptibilities, 
as an infinite sum of D-finite functions, should be expected at first sight. 
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for all their singularities are rational numbers, and their Wronskians are iV-th roots 
of rational functions [24] . Furthermore, it has been shown later that these x^ ' s are ; 
in fact, solutions of globally nil-potent operators [25], or G-operators [26, 27]. It is 
worth noting that the series expansions, at the origin, of the X Sj m a well-suited 
variable [13, 23] w, actually have integer coefficients, even if this result does not have 
an immediate prooff for all integers n (in contrast with the full susceptibility). From 
the first truncated series expansions of x^ n \ the coefficients for generic n can be 
inferred [25] 

X (n) (w) = 2"-w™ 2 ■ (l + An 2 -w 2 + 2-(4n 4 +13n 2 + 1) ■ w 4 

+ ^.(n 2 +4)(4n 4 + 23n 2 + 3)-w 6 (1) 
+ i • (32 n 8 + 624 n 6 + 4006 n 4 + 8643 n 2 + 1404) • w 8 

O 

+ — • (n 2 + 8)- (32 n 8 + 784n 6 + 6238?i 4 + 16271 n 2 + 3180) • w 10 +■■■). 
15 / 

Note that these coefficients of (w) / 2 n , which depend on n 2 , are integer coefficients 
when n is any integer, this integrality property of the coefficients for any integer n 
being not straightforward (see [17]). These coefficients arc valid up to w 2 for n > 3, 
w for n > 5, w e for n > 7, w 8 for n > 9, and w 10 for n > 11 (in particular 
it should be noted that x lS an even function of w only for even n). Further 
studies on these x^'s showed the fundamental role played by the theory of elliptic 
functions^] (elliptic integrals, modular forms) and, much more unexpectedly, Calabi- 
Yau ODEs [28, 29]. These recent structure results thus suggest to see the occurrence 
of series with integer coefficients as a consequence of modularity [30] (modular forms, 
mirror maps [28, 29, 30, 31], etc) in the Ising model. 

Along this line, many other examples of series with integer coefficients emerged 
in mathematical physics (differential geometry, lattice statistical physics, enumerative 
combinatorics, replicable functions^ ...). One must, of course, also recall Apery's 
results [37]. Appendix A of [17] gives a list of modular forms, and their associated 
series with integer coefficients, corresponding to various lattice Green functions [38, 39, 
40, 41], that are, often, expressed in terms of HeunG functionsf which can be written 
as hypcrgcomctric functions with two alternative pullbacks [17] (see also sections (6.1) 
and (6.2) below). Let us underline, in Appendix A of [17], the Green function for the 
diamond lattice [41], the Green function for the face- centred cubic lattice (see equation 
(19) in [41]), and more examples corresponding to the spanning tree generating 
functions [42] (and Mahler measures). This integrality is also seen in the nome and 
in other quantities like the Yukawa coupling [28]. 

In this paper we restrict on series with integer coefficients, or, more generally, 
globally bounded [43] series of one complex variable, but it is clear that this integrality 
property does also occur in physics with several complex variables: they can, for 

t We are interested in this paper in the emergence of integers as coefficients of D-finite series. In 
general, this emergence is not obvious: it cannot be simply explained at the level of the linear 
recurrence satisfied by the coefficients, as illustrated by the case of Apery's calculations (see also 
section (6.2) and Appendix G in [17]). 

If Which is not a surprise for Yang-Baxter integrability specialists. 

(j The concept of replicable functions is closely related to modular functions [32] , (see the replicability 
of Hauptmoduls). Calabi-Yau threefolds, and more generally the concept of modularity [30, 33, 34, 
35, 36]. 

f Gcnerically HcunG functions are far from being modular forms. 
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instance, be seen for the previous (D-finite§) n-fold integrals x f° r the anisotropic 
Ising model [44] (or for the Ising model on the checkerboard lattice) , or on the example 
of the lattice Ising models with a magnetic field! (see for instance, Bessis et al. [4]). 

One purpose of this paper is to "disentangle" the notion of series with integer 
coefficients (integrality) and the notion of modularity [30, 33, 34, 35, 36]. In this "down- 
to-earth" paper we will use the wording of "modularity" , not to refer to the modularity 
conjecture and other Serre's results that certain Geometric Galois representations are 
modular, but as a quick proxy word to say that a series solution of a linear differential 
operator as well as the nome, and hopefully other series (Yukawa coupling, ...) are all 
series with integer coefficients. 

We will show that the x s are series globally bounded series, as a consequence 
of the fact that they actually are diagonals of rational functions for any value of 
the integer n. We will generalise these ideas, and show that an extremely large 
class of problems of mathematical physics can be interpreted in terms of diagonals of 
rational functions: n-fold integrals with algebraic integrand of a certain type that we 
will characterise, Calabi-Yau ODEs, MUM linear ODEs [46], scries whose coefficients 
are nested sums of binomials, etc. We take, here (see also [17]) a learn- by-examplc 
approach: on such questions one gets a much deeper understanding from highly non- 
trivial examples than from general mathematical demonstrations [47]. 

2. Series integrality, diagonal of rational functions 

Let us recall some concepts that will be fundamental in this paper, first the notion 
of globally bounded series, and, then, the concept of diagonal of a "function"^], and 
some of its most important properties. The main reason to introduce this concept 
of diagonal of function, not very familiar to physicists, is that it enables to consider 
diagonal of rational functions, this class of functions filling the gap between algebraic 
functions and G-scrics: they can be seen as generalisations of algebraic functions. 
Thus this class of functions can play a key role to decipher the complexity of functions 
occurring in theoretical physics. 

2.1. Globally bounded series 

Let us first recall the definition of being globally bounded [43] for a series. Consider a 
series expansion with rational coefficients, with non-zero radius of convergencef . The 
series is said to be globally bounded if there exists an integer N such that the series 
can be recast into a scries with integer coefficients with just one rescaling x — > N x. 

A necessary condition for being globally bounded is that only a finite number of 
primes occur as factors of the denominators of the rational number series coefficients. 
There is also a condition on the growth of these denominators, that must be bounded 
exponentially [43], in such a way that the series has a non-zero p-adic radius of 
convergence for all primes p. When this is the case, it is easy to see that these 
series can be recast, with just one rescaling, into series with integer coefficientsi. 

§ For several complex variables the ODEs of the paper are replaced by Picard-Fuchs systems. 

J Along this line original alternative representations of the partition function of the Ising model in a 

magnetic field are also worth recalling [45]. 

If This is an abuse of language: the "functions" are in fact defined by series of several complex 

variables: they have to be Taylor, or Laurent, series (no Puiseux series). 

f A series like the Euler-series JZ^Lo n ' ' x " wruc h has integer coefficients is excluded. 

J For a first set of series with integer coefficients, see appendix A of [17], where a set of such series 
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2.2. Definition of the diagonal of a rational function 

Assume that T(z\, . . . , z n ) = P(z\, . . . , z n )/Q(zi, . . . , z n ) is a rational function, 
where P and Q are polynomials of Z\ , • • • , z„ with rational coefficients such that 
Q(0, . . . , 0) 7^ 0. This assumption implies that T can be expanded at the origin as a 
Taylor series with rational number coefficients 

oo oo 

T(z U Z 2 , Z„) = E E F mi,..,m n - Z? 1 ■■■ Z%». (2) 

mi — m n — 

The diagonal of T is defined as the series of one variable 

oo 

DiagLFLzi, z 2 , 3 n )) = E Fm >™> ->™ ' z ™- ( 3 ) 

m = 

More generally, one can define, in a similar way, the diagonal of any multivariate 
power series J 7 , with rational number coefficients, or with coefficients in a finite field|. 



2. 3. Main properties of diagonals 

The concept of diagonal of a function has a lot of interesting properties (sec for 
instance [49]). Let us recall, through examples, some of the most important ones. 

The study of diagonals goes back, at least, to Polya [50], in a combinatorial 
context, and to Cameron and Martin [51] in an analytical context related to Hadamard 
products [52] . Polya showed that the diagonal of a rational function in two variables 
is always an algebraic function. The most basic example is T = 1/(1 — Z\ — 22), for 
which 



Diag(^) = Diag £ £ 

\mi — ni2— 



m 1 + m 2 
mi 



m— v 7 v 

The proof of Polya's result is based on the simple observation that the diagonal 
Diag(J r ) is equal to the coefficient of z\ in the expansion of Tizx^z f z\). Therefore, 
by Cauchy's integral theorem, Diag(J r ) is given by the contour integral 

Diag(^) = -L I^ ZliZ / Zl )^l ) ( 5 ) 
2m J 1 z\ 

where the contour 7 is a small circle around the origin. Therefore, by Cauchy's 
residue theorem, Diag(J r ) is the sum of the residues of the rational function Q = 
T(z\ y zj z\)l Z\ at all its singularities s(z) with zero limit at z = 0. Since the residues 
of a rational function of two variables are algebraic functions, Diag(J r ) is itself an 
algebraic function. 

For instance, when J- = 1/(1 — z\ — Z2), then Q = T(z\,zj z\)j Z\ has two 
poles at s = 5(1 ± y/1 — 4z). The only one approaching zero when z — s- is 
s = |(1 — \/l — 4z). If p(s)/q(s) has a simple pole at So, then its residue at so is 
p(so) /q' (so) . Therefore 

Diag(^) = J- /— ^ = = -j^. (6) 

' 2m J zi-zf-z 1 - 2s VT^Tz V ; 



with integer coefficients corresponding to modular forms is displayed. See also (6.1) and (6.2) below. 
\ The definition even extends to multivariate Laurent power series, see e.g. [48] . 
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2.3.1. Diagonals of rational functions of more than two variables 

When passing from two to more variables, diagonalisation may still be interpreted 
using contour integration of a multiple complex integral over a so-called vanishing 
cycle [53]. However, the result is not an algebraic function anymore. A simple example 
is J- = 1/(1 — Z2 — z-t, — z\Z2 — Z1Z3), for which 

Diag( J") = 1 + 4z + 36z 2 + 400z 3 + 4900z 4 + 63504z 5 + • • • (7) 
is equal to the complete elliptic integral of the first kind 



Diag(J-) = Err) = 2 F 1 ([i,i],[l];16z), (8) 

m>0 ^ ' 

which is a transcendental function. A less obvious example (see [54] for a related 
example with a combinatorial flavor) is 

Diag (- ) (9) 

V 1 - Zi - z-z - z 3 - ZiZ 2 - z 2 z 3 - z 3 zi - z 1 z 2 z 3 J 

1 (I 2 54z 

r-^' 2Fl l [ 3'3 ul]; (1^7)3, 

It was shown by Christol [55, 56, 57] that the diagonal Diag(J r ) of any rational 
function T is D-finite, in the sense that it satisfies a linear differential equation 
with polynomial coefficients^]. Moreover, the diagonal of any algebraic power series 
with rational coefficients is a G-function coming from geometry, i.e. it satisfies the 
Picard-Fuchs type differential equation associated with some one-parameter family 
of algebraic varieties. Diagonals of algebraic power series thus appear to be a 
distinguished class of G-functionsjj. It will be seen below (see (2.5)) that algebraic 
functions with n variables can be seen as diagonals of rational functions with 2 n 
variables. Thus diagonals of rational functions also appear to be a distinguished class of 
G-functions. It is worth noting that this distinguished class is stable by the Hadamard 
product: the Hadamard product of two diagonals of rational functions is the diagonal 
of rational function. 

An immediate, but important property of diagonals of rational functions, with 
rational number coefficients, is that they are globally bounded, which means that they 
have integer coefficients up to a simple change of variable z — > N z, where N G Z. 

2.3.2. Diagonals of rational functions modulo primes 

Furstenberg [60] showed that the diagonal of any multivariate rational power series 
with coefficients in a field of positive characteristic is algebraic. Dcligne [53, 48] 
extended this result to diagonals of algebraic functions. For instance, when J- = 
1/(1 — Z2 — Z3 — Z1Z2 — Z1Z3), one gets modulo 7 

Diag(J") mod 7 = 1 + 4 z + z 2 + z 3 + Az 7 + 2z 8 + 4z 9 +••• 

= 1 mod 7. (10) 

#1 +4z +z 2 +z 3 

More generally, in this example, for any prime p, one has 

Diag(J") = P(z) 1/(1 - p) modp (11) 

A more general result was proved by Lipshitz [58] : the diagonal of any D-finite series is D-finite, 
see also [59]. 

tt Such diagonals are solutions of G-operators. They are functions that are always algebraic mod. 
a prime p. They fill the gap between algebraic functions and G-series: they can be seen as 
generalisations of algebraic functions. 
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where the polynomial P{z) is nothing, but [61, 62, 63] 

P(*)= 2^1 ([^], [1];16^) 1_P modp= £ ( P ^ /2 ) -( 16z )"- ( 12 ) 

n=0 ^ ' 

Note, however, that the Furstensberg-Dcligne result [60, 53], that we illustrate, 
here, with J- = 1/(1 — z 2 — z% — z\z 2 —Z1Z3), goes far beyond the case of hypergeometric 
functions for which simple closed formulae can be displayed. 

2. 4- Hadamard product and other products of series 

Let us also recall the notion of Hadamard product [52, 64] of two series, that we will 
denote by a star. 



oo 



If f(x) = ^ «„ ■ i", g(x) = ^2b n -x n , then: 

n=0 n=0 
oo 

f(x)*g(x) = ^On-bn-X*. (13) 

n=0 

The notion of diagonal of a function and the notion of Hadamard product are 
obviously related: 

Diagf /i(x x ) • f 2 {x 2 ) ■ ■ ■ /«(£«)) = fi(x) * f 2 {x) * ■ ■ ■ * f n (x). (14) 

In other words, the diagonal of a product of functions with separate variables is equal 
to the Hadamard product of these functions in a common variable. In particular, 
the Hadamard product of n rational (or algebraic, or even D-fmitc) power series is 
D-fmite§. 

The Hadamard product of two scries with integer coefficients is straightforwardly 
a series with integer coefficients. Furthermore, the Hadamard product of two operators, 
annihilating two series, defined as the (minimal order, monic) linear differential 
operator annihilating the Hadamard product of these two series, is a product compatible 
with a large number of structures and concepts^ that naturally occur in lattice 
statistical mechanics. We have a similar compatibility property between the diagonal 
and the Hurwitz product [17, 65]. 

2.5. Furstenberg's result on algebraic functions 

It was shown by Furstcnbcrg [60] that any algebraic series in one variable can 
be written as the diagonal of a rational function of two variables. The basis of 
Furstenberg's result is the fact that if f(x) is a power series without constant term, 
and is a root of a polynomial P(x,y) such that P y (0, 0) ^ 0, then 

'M " Diag {" 2 ' pisri/l) whe " : p > ^ f' <15) 

When P y (0,0) = 0, formula (15) is not true anymore. However, Furstenberg's result 
still holds [17]. 

§ The Hadamard product of rational power series is still rational, but the Hadamard product of 
algebraic series is in general transcendental. 

| For instance, the Hadamard product of two globally nilpotent [25] operators is also globally 
nilpotent. 
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Note that this representation as diagonal of a rational function is, by no means 
unique, as can be seen on the algebraic functionf 

X 1 n 3 q O A 35 e 63 c . . 

f= vr^ = x + r + i x + t, x +Y2* x +w, x + (16) 

which is the diagonal of (2xy — cx + cy)/(x +y + 2) for any rational number c. 

Furstcnbcrg's proof does not necessarily produce the simplest rational function 
(sec [17]). 

Furstcnbcrg's result has been generalised to power series expansions of algebraic 
functions in an arbitrary number of variables n: any algebraic power series^ with 
rational coefficients is the diagonal of a rational function with 2n variables (see Denef 
and Lipshitz [66]). 

3. Selected n-fold integrals are diagonals of rational functions 

Among many multiple integrals that are important in various domains of mathematical 
physics, and before considering other n-fold integrals of the "Ising class" , let us first 
consider the n-particle contribution to the magnetic susceptibility of the Ising model 
which we denote x ( w )- They are given by (n — l)-dimensional integrals [11, 67]: 

*«<-> - ^(nft)--^i <«> 

j=l J o 

where, defining <t> by J2i=o ®i = 0> we se t 
G = FT 2 - »«»(♦«-»<) where . A . = j _2„cos(4 ( ). (18) 

0<i<j<n-l y 101 

The integrality property (1) had been checked [12] for the first x^'s and inferred [25] 
for generic n. We are going to prove it\ for any integer n, showing a much fundamental 
result, namely that all the {n—Vj-fold integrals 's are very special: they are actually 
diagonals of rational functions. 

3.1. x^ as a t°y example 

At first sight the x s are involved transcendental holonomic functions. Could it be 
possible that they correspond to the distinguished class [48] of G-functions, generalising 
algebraic functions, which have an interpretation as diagonals of multivariate algebraic 
functions (and consequently diagonals of rational functions with twice more variables)? 
If this is the case, then the series of the x a mu st necessarily reduce modulo any 
prime to an algebraic function (see (2.3.2)). The x^ an d contributions being 
too degenerate (a rational function and a too simple elliptic function), let us consider 

f Here, / is annihilated by P(x, y) = (1 — x)y 2 — x 2 , which is precisely such that P y {0, 0) = 0. 
K In the one-variable case, Puiseux series could be considered but only after ramifying the variable, 
f For the x' n 'i the rescaling factor ("Eisenstein constant") is 2 or 4 according to the fact that one 
considers high or low temperature series [12, 24]. 
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the first non-trivial case, namely x ■ Its series expansion has already been displayed 
in [11]. It reads x (3) / 8 = ^ 9 ' F ( w ) with: 
F(w) = 1 + 36 w 2 + 4w 3 + 884 w 13 + 196 w 5 + 18532 w 6 + 6084 w 7 +••• 

Since we have obtained the exact ODE satisfied by x^ we can produce as many 
coefficients as we want in its series expansion. Let us consider this series modulo the 
prime p = 2. It now reads the lacunary series 

F(w) mod 2 = 1 + w 8 + w 24 + w 56 + w 120 + w 248 + w rM4 + w 1016 + • • • , 

solution of the functional equations on F(w) or, with z = w 8 , on G(z) = 
I + w s ■ F(w) 

F(w) = 1 + w 8 ■ F(w 2 ), G{z) = z + G(z 2 ), (19) 

where one recognises, with equation G(z) = z + G(z 2 ), Furstcnberg's example [60] 
of the simplest algebraic function in characteristic 2§. In fact H(w) = w 9 F(w) is 
solution of the quadratic equation: 

H{w) 2 + w ■ H(w) +w 10 = mod 2. (20) 

The calculations are more involved modulo p = 3. Indeed, H(w) = x^( w )/^ 
satisfies, modulo 3, the polynomial equation of degree nine 

p 9 ■ H(w) 9 + w 6 ■ p 3 ■ H{wf + w 10 ■ pi ■ H(w) 

+ w 19 ■ • p^ = 0, (21) 

where: 

Pg = (w + l) 3 (w 2 + l) 18 (w - l) 24 , 

p 3 = {w 2 + l) 18 (1 - w) 15 (w 4 - w 2 - 1), p x = (w 2 + l) 20 (l-«;) 13 , 

P { q ] = w 6 + w 5 +w 4 -w 2 -W + l, (22) 

p { V = w 37 - w 36 + w 35 ~ w 33 + w 31 - w 30 + w 28 + w 27 + w 24 - w 23 + w 22 

- w 21 - w 18 - w 16 + w 14 - w 12 - w 11 - w w + w 7 - w 5 - w 3 - 1. 

The calculations are even more involved modulo larger primes. The (minimal 
order) linear differential operator annihilating the series mod. 5, readsf: 

(x+l)(a; 2 + a; + l)(a; + 2)- x 4 ■ D 4 +2x 3 - {x 3 + 2 x 2 + Ax + 4) (x + 4) • D 3 X 

+ x 2 ■ (x 4 + 3 x 3 + 4) • D 2 X + 4 • {x 4 + 3) • x ■ D x +3 (23) 

If one can easily get this linear differential operator, finding the minimal polynomial 
of modulo 5, generalising (20) or (21), or rather, the polynomial P{k, w), where 
K = x^ 3 ^(w)/w 9 , such that P(«, w) = mod. 5, requires a very large number of 
coefficients. The polynomial [17] P(k, w) is of degree 50 in k, of degree 832 in x and 
is the sum of 4058 monomials. 

One can imagine, in a first step, that the series mod. any prime p are also 
algebraic functions, and, in a second step, that x^ m ay be the diagonal of a rational 
function. In fact we are going to show, in the next section, a stronger result: the 
X^™' 's are actually diagonals of rational functions, for any integer n. 

§ Modulo the prime p = 2, the previous functional equation becomes G(z) = z + G(z) 2 . 
f This operator is of zero 5-curvature [25]. 
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3.2. The x 's are diagonals of rational functions 

Let us, now, consider the general case where n is an arbitrary integer. With the 
change of variable z k = exp(i$fe) (where i 2 = —1), one clearly gets 

I[zk = 1, = id*j , (24) 

fe=o Z] 

2cos($ fc ) = z k +— , 2cos($ fc -$ i ) = ^ + ^-, 

Zjz Zj Zk 

and (17) becomes 

- ^(n^if)-^.- t») 

where C is the path "turning once counterclockwise around the unit circle" and where 
F is algebraic over Q(w, Z\, . . . , z n _i) and reads: 

F(w, Zl ,...,z n -x) = Y ■ X"" 1 • • G. (26) 

Now, let us suppose that F is analytic] at the origin, namely that it has a Taylor 
expansion (2). Then applying (n — 1) times the residue formula, one finds 

X (n) H = Diag( (2zo -^"- l) " • F(z ••• Zn-x, *i, .... *,_!)). (27) 

To check that this is actually true, we introduce an auxiliary set, namely T n the subset 
of Laurent series Q[zi, . . . , 1 , . . . , z^-J^w]], consisting of series 

oo 

/(«>, Zi, . . . , Z n -x) = P m - w m , 

m=0 

where P m belongs to Q[z\, . . . , z n -i, z^ 1 , . . . , z~\] and is such that the degree of 
P rn , in each of the z^ 1 , is at most m. 

Then to prove that F(zqz\ . . . z n —i, Zi, . . . , z n -i) has a Taylor expansion, we only 
have to verify that F(w, z\, . . . , z n -i) belongs to this auxiliary set T n . Checking this 
is a straightforward stcp-by-stcp computation on auxiliary functions: 

A k = 1 -w (z k +^)> for: l<i<n-l, 

A k = 1 - w ■ ( h z\- ■ ■ Zn-i) , for: i = 0. 

\zi ■ ■ ■ z„_i J 

Hence A k belongs to this auxiliary set T n . So to be sure that the inverse or the 
square root of some function in this auxiliary set T n is also in this auxiliary set T n 
we only have to check that its first Taylor coefficient is actually 1, or w, w 2 , or w n , 
n integer. It is straightforward to see that: 

A 2 -Aw 2 = 1 -2w(z k +—) + w 2 .(z k -—) , (28) 



hence: J A\ - Aw 2 = 1 + • • • , 

f One could consider Laurent, instead of Taylor, expansions, but this is a slight generalisation [48, 
68, 69]. 
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Vk = * = „ = 1 + • ■ • , Y = 1 + 

2w 

Xh = = W + • • ■ , Xl. X; = W 

A k + y/Al - 4^ 

X = w n + • • • , — — = !+•••, 



2 



G* 



-(Z fc -Zj) 2 \ T T f-(z k -Zj) 



n (:1 .r...r,; 2 ••,.:.) ...II 



. (1 Xf~Xj)^ • zj~ Zj J V z/^ Zy 



Thus, this shows that F belongs to the auxiliary set T n - 

F(w,zi,...,z n -i) = (29) 

oo Mi(m) M„_i(m) 

E\ " \ " - „,m .mi+m m„_i+m 

/ y '' / y «m,m 1 ,...,m n _i <" ■■• ^n-1 

m—0mi — — m m n -i— — m 

Consequently, it makes sense to take its diagonal. The residue theorem requires 
searching the terms not containing Zj i.e. such as mi = . . . = m n _i = 0. One 
therefore gets: 

oo 

x {n) = £ a ^o,...,o • w ™- ( 3 °) 

m=0 

In particular, we find 

oo Mi(m) M„_i(m) 

x (n) = Dia g (]T £ •■• E w,™,-!- ^r +ra ---c , r 1+m 

m— mi--m m^-i— — m 

= Diag(V(z zi . . . z„_i, zi, . . . , (31) 

As F is algebraic, jj; is the diagonal of an algebraic function of n variables 
and, consequently, the diagonal of a rational function of 2 n variables. 

We thus see that we can actually find explicitly the algebraic function such 
that its diagonal is the n-fold integrals x ■ it is nothing but the integrand of 
the n-fold integral, up to simple transformations (namely F(w, zi, z„_i) — > 
F(z zi ... z„_i, z x , z„_i)). 

Remark : x'™' is a solution of a linear differential equation, and has a radius of 
convergence equal to 1/4 in w. Among the other solutions of this equation, there is 
the function obtained by changing the square root appearing in x k into its opposite. 
A priori there are 2 n ways to do this, hence 2™ new solutions but, not all distinct. At 
first sight, for these new solutions, the Xfc's are no longer in 7~ n . 

In fact, we find some quite interesting structure. Let us consider, for instance, the 
case of x^ • If one considers other choices of sign in front of the nested square roots 
in the integrand, the series expansions of the corresponding n-fold integrals read: 

w + 6 w 2 + 28 w 3 + 124 + 536 w 5 + 2280 w 6 + 9604 w 7 + 40164 w 8 

+167066 w 9 + 692060 w 10 + 2857148 w 11 +■■■ 
w 2 + 6 w 3 + 30 w 4 + 140 w 5 + 628 w 6 + 2754 w 7 + 11890 w 8 + 50765 w 9 

+214958 w 10 + 904286 w 11 +■■■ 
One does remark that all these alternative series are, as x' 3 \ series with integer 



coefficients [17]. 
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3.3. More n-fold integrals of the Ising class and a simple integral of the Ising class 

It is clear that the demonstration we have performed on the X s can a l so be 
performed straightforwardly, mutatis mutandis, with other n-fold integrals of the 
"Ising class" like the n-fold integrals <&h in [67] , which amounts to getting rid of the 
fermionic term G (see (26)), the Xd"''' s corresponding to n-fold integrals associated 
with the diagonalf susceptibility [15, 29] (the magnetic field is located on a diagonal 
of the square lattice), the 3>£^'s in [14] which are simple integrals, and also for all the 
lattice Green functions displayed in [41, 46], and the list is far from being exhaustive. 
For instance, the simple integral is the diagonal of the algebraic function: 



4 • (1 - t 2 y 1/2 ^P - A , where: (32) 

F n = 1 -2w +(l-4w + Aw 2 -4w 2 t 2 ) 1/2 , (33) 

G n = 1 -2wt-T n _ 1 (^j +((l -2wt- T n -i(j)) 2 -4w 2 ■ t 2 y /2 , 

and where T„_i(i) is the (n — l)-th Chebyshev polynomial of the first kind. The way 
we have obtained these Chebyshev results (32) is displayed in [17]. 

As opposed to the x^'s, the integral (iy) * s the diagonal of an algebraic 
function of two variables (and, thus, the diagonal of a rational function of four 
variables) independently of the actual value of n. 



A\x, zi, z 2 , ■■■ , (34) 



3.4- More general n-fold integrals as diagonals 

More generally the demonstration we have performed on the X s can be performed 
for any n-fold integral that can be recast in the following form: 

f f f dzi dzi dz n 

Jc JO JO z l z 2 Z n 

where the subscript C denotes the unit circle, and where A denotes an algebraic 
function of the n variables, which (this is the crucial ingredient), as a function of 
several variables x and the Zfc's, has an analytical expansion at (x, Z\, z%, ■ ■ ■ , z n ) = 
(0, 0, 0, ••• , 0): 

A\x, zi, Z2, ■■■ , Zn^j = (35) 

oo oo oo oo 

E\ " \ " . . . \ " A . ,mi m 2 _ _ m„ . m 

/ j / j / t n m,mi,m2,-,m„ 6\ *2 n A ■ 

m — mi — m2 — m n — 

Consequently, an extremely large set of n-fold integrals occurring in theoretical 
physics (lattice statistical mechanics, enumerative combinatorics, number theory, 
differential geometry, ...) can actually be seen to be diagonals of rational functions. 
These n-fold integrals correspond to series expansions (in the variation parameter x) 
that are globally bounded (can be written after one rescaling into series with integer 
coefficients), and are solutions of globally nilpotent [25] linear differential operators. 

f Of course this "diagonal [15, 29] wording" should not be confused with the notion of diagonal of a 
function. 
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Such a general n-fold integral is, thus, the diagonal of an algebraic function (or 
of a rational function with twice more variables [66]) which is essentially the integrand 
of such n-fold integral. Furthermore, such a general n-fold integral is solution of a 
(globally nilpotent [25]) linear differential operator, that can be obtained exactly from 
the integrand, using the creative telescoping method [17]. 

Finally, in the case of Calabi-Yau ODEs (see below), these functions can 
be interpreted as periods of Calabi-Yau varieties, these algebraic varieties being 
essentially the integrand of such n-fold integrals. The integrand is thus the key 
ingredient to wrap, in the same bag, the algebraic geometry viewpoint, the differential 
geometry viewpoint and analytic and arithmetic approaches (series with integer 
coefficients). 



4. Calabi-Yau ODEs solutions and series with binomials seen as diagonals 

Recalling Calabi-Yau ODEs 

Calabi-Yau ODEs have been defined in [70] as order-four linear differential ODEs 
that satisfy the following conditions: they are maximal unipotent monodromy [71, 72] 
(MUM), they satisfy a "Calabi-Yau condition" which amounts to imposing that the 
exterior squares of these order-four operators are of order five (instead of the order 
six one expects in the generic case), the series solution, analytic at x = 0, is globally 
bounded (can be reduced to integer coefficients), the series of their nome and Yukawa 
coupling arc globally bounded. In the literature, one finds also a cyclotomic condition 
on the monodromy at the point at oo, x = oo, and/or the conifoldf character of one 
of the singularities [74] . 

Let us recall that a linear ODE has MUM (maximal unipotent monodromy [28, 73]) 
if all the exponents at (for instance) x = are zero. In a hypergeometric framework 
the MUM condition amounts to restricting to hypergeometric functions of the type 
n+iF n ([a±, 02, • • • On]i [1, 1, ■ ■ • 1], x), since the indicial exponents at x = are the 
solutions of p (p + b\ — 1) ■ • • (p + b n — 1) = p n+1 = 0, where the bj are the lower 
parameters which are here all equal to 1. 

Let us consider a MUM order-four linear differential operator. The four solutions 
yo, 2/i, 2/2, 2/3 of this order- four linear differential operator read: 

ln(a;) 2 

yo, yi = yo- M 35 ) + h, yi = vo — = h ^ ' ln ( x ) + 

ln(a;) 3 ln(x) 2 

2/3 = yo — g — + yi — 2 — + 2/2 ' ( x > + y3 ' 

where j/o, 2/1, 2/2, 2/3 are analytical at x — (with also yi(0) = 2/2(0) = 2/3(0) = 0). 
The nome of this linear differential operator reads: 

q(x) = exp( — ] = x- expf— V (36) 
^ 2/0 ' \yo> 

Calabi-Yau ODEs have been defined as being MUM, thus having one solution 
analytical at x = 0. As far as Calabi-Yau ODEs are concerned, the fact that this 
solution analytical at x — has an integral representation, and, furthermore, an 
integral representation of the form (34) together with (35), is far from clear, even if 



f The local exponents are 0,1,1,2. For the cyclotomic condition on the monodromy at 00, see 
Proposition 3 in [73]. 
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one may have a "Geometry-prejudice" that this solution, analytical at x = 0, can be 
interpreted as a "Period" and "Derived From Geometry" [26, 27, 75]. 

Large tables of Calabi-Yau ODEs have been obtained by Almkvist et al. [73, 
76, 77]. It is worth noting that the coefficients A n of the series corresponding to 
the solution analytical at x = 0, arc, most of the time, nested sums of product 
of binomials, less frequently nested sums of product of binomials and of harmonic 
numbers^ H n , and, in rare cases, no "closed formula" is known for these coefficients. 

Let us show, in the case of A n coefficients being nested sums of product of 
binomials, that the solution of the Calabi-Yau ODE, analytical at x = 0, which 
is by construction a series with integer coefficients, is actually a diagonal of rational 
function, and furthermore, that this rational function can actually be easily built. 



4-2. Calculating the rational function for nested product of binomials 

For pedagogical reasons we will just consider, here, a very simple example§ of a series 
S(x), with integer coefficients, given by a sum of product of binomials 

OO Tl / \ 3 

5 ( s ) = EEC = HeunG(-l/8,l/4,lAAA;-x) (37) 

n=0 k=0 ^ ' 

= 1 +2x +10 a; 2 + 56 x 3 + 346 x 4 + 2252 a; 5 + 15184 a; 6 + ■■■ 

This is the generating function of sequence A in Zagier's tables of binomial coefficients 
sums (see p. 354 in [78]). 

The calculations of this section can straightforwardly (sometimes tediously) be 
generalised to more complicated [79] nested sums of product of binomialsf. 

Finding that a series is a diagonal of a rational function amounts to framing it 
into a residue form like (34). In order to achieve this, we write the binomial (^) as 
the residue 

^ ■ "A (38) 

z K z 

and, thus, we can rewrite S(x) as 
(2i7r) 3 ■ S(x) = 

(n+1) 

, n dz\ dz2 dz 3 







i 




G) 




2iir 


L 



OO 

E- 

n=0 



- [l/{ziz 2 z 3 )\ , 
(l+z 1 )(l + z 2 )(l + ^ 3 )- x) 

l-(l/(*l*2*3)) V ' ZlZ2Z3 

V' z i z 2 Z3 (,-. W1 . Wl . x \" dz 1 dz 2 dz 3 

>- I (1 + Zl) (1 + Z 2 ) (1 + Z 3 ) ■ X 1 • 

1 - Zi Z 2 Z 3 \ J ZiZ 2 Z 3 

n—0 
oo 1 

y — - 

^ 1 - Zl 



(1 + zx) (1 + z 2 ) (1 + z 3 ) ■ x\ n dzi dz 2 dz 3 



zi z 2 z 3 \ Zl z 2 z 3 / zi z 2 z 3 

% The generating function of Harmonic numbers is H(x) = H n ■ x n — — ln(l — x)/(l — x). 
§ See Proposition 7.3.2 in [72]. 

f Not necessarily corresponding to modular forms as can be seen on (47), (48). 
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ul s d Zl dz 2 dz 3 

R(x; zi, z 2 , z 3 ) ■ , (39) 

Z\ z 2 z 3 

where R(x; z±, z 2 , z 3 ) reads: 

zi z 2 z 3 

(1 - x- (1 + Z 1 )(1 + 2 2 )(1 + Z 3 )) (z lZ2 z 3 - x- (1 + Zl )(l + z 2 )(l + z 3 ))' 

From this last result one deduces immediately that (37) is actually the diagonal of: 

1 

(1 - 2 • (1 + zi)(l + z 2 )(l + z 3 )) • (1 - z zi z 2 z 3 (1 + Zl )(l + z 2 )(l + 2:3)) ' 

Note that, as a consequence of a combinatorial identity due to Strchl and 
Schmidt [80, 81, 82], S(x) can also be written as 

*) - t 1(:) 2 ( 2 ;)-" - 1 1 G)'c,V- <4o) 

n=0 fe=0 v 7 v 7 n=0 fe=[n/2] v 7 v 7 

Calculations similar to (39) on this alternative binomial representation (40), enable 
to express (37) as the diagonal of an alternative rational function: 

(41) 

(1 -zq> (l + z 1 )(l + z 2 )(l + z 3 ) 2 ) • (1 -z oZl z 2 - {l + Zl )(l + z 2 ))' 

We thus see that we can actually get explicitly, from straightforward calculations, 
the rational function (40) for the Calabi-Yau-like ODEs (occurring from differential 
geometry or enumerative combinatorics) when series with nested sums of binomials 
take place, and, more generally, for enumerative combinatorics problems (related or 
not to Calabi-Yau manifolds) where series with nested sums of binomials take place. 

These effective calculations are actually algorithmic, and guarantee to obtain an 
explicit expression for the rational function (40). However the rational function is far 
from being unique, and worse, the number of variables, the rational function depends 
on, is far from being the smallest possible number. Finding the "minimal" rational 
function (whatever the meaning of "minimal" may be) is a very difficult problem. 
Appendix A provides a non-trivial illustration of this fact with explicit calculations 
on the well-known Apery series and its rewriting due to Strehl and Schmidt [80, 81, 82]. 
We see in a crystal clear way in Appendix A that, when a given function is a diagonal 
of a rational function, the rational function is far from being unique, the "simplest" 
representation (minimal number of variables, lowest degree polynomials, ...) being 
hard to find. Similar computations show that the generating function of sequence B 
and E in Zagier's list [78] are both diagonals of rationals function in four variables. 
These results are given in section (5.1) of [17]. 

All these calculations can systematically be performed on any series defined by 
nested sums of product of binomials. We have performed such calculations on a large 
number of the series corresponding to the list of Almkvist et al [73] , that are given by 
such nested sums of product of binomials. 

5. Comments and speculations 

5.1. Christol's theorem 



In [43] (page 61 Theorem 12, see also Proposition 7 in page 50 of [57]) it is proved 
that any power series with an integral representation and of maximal weight for the 
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corresponding Picard-Fuchs linear differential equation is the diagonal of a rational 
function and, in particular, is globally bounded. 

The technical nature of the original papers is such that the result itself is difficult 
to find. This paragraph is devoted to explain, in down-to-earth terms, the somewhat 
esoteric expressions used in its wording, and to explain what it means on explicit 
examples. As the original proof is very obfuscated its principle is sketched in [17]. 

Disappointingly, when applied to a hypcrgeometric n +iF n , this result becomes 
somewhat trivial. More precisely, the hypcrgeometric function is of maximal weight 
if and only if bj = 1 for all j (there is only n!'s in the denominator of coefficients). 
In that case it is obviously the Hadamard product of algebraic functions, therefore 
diagonal of a rational function: 

n F n -i([ai, a 2 , ■ ■ ■ , On], [1, 1, • • • 1], x) (42) 
= (1 - x y ai * (1 - x)' a2 ■■■ * (1 - a;) - "". 

Therefore, we now have (at least) three sets of problems yielding diagonal of 
rational functions: the n-fold integrals of the form (34) with (35), the Picard-Fuchs 
linear ODEs with solution of maximal monodromy weight and, finally, the problems 
of enumerative combinatorics where nested sums of products of binomials take place. 
Diagonal of rational functions, thus, occur in a quite large set of problems of theoretical 
physics. 

5.2. Christol's conjecture 

The diagonal of a rational function is globally bounded (i.e. it has non zero radius of 
convergence and integer coefficients up to one rescaling) and D-finite (i.e. solution of 
a linear differential equation with polynomial coefficients)^:. 

The converse statement is "Christol's conjecture" [43] saying that any D-finite, 
globally bounded series is necessarily the diagonal of a rational function. 

A fantastic Chudnovski theorem ([84] page 267) asserts that the minimal 
linear differential operator of a G-function (and in particular of a D-finite globally 
bounded series) is a G-operator (i.e. at least conjecturally, a globally nilpotent 
operator) [25, 26, 27]. "Christol's conjecture" amounts to saying something more: 
if the solution of this globally nilpotent linear differential operator is, not only a 
G-series, but a globally bounded series, then it is the diagonal of a rational function. 

Conversely the solution, analytical at 0, of a globally nilpotent linear differential 
operator is necessarily a G-function [26, 27]. Moreover, a "classical" conjecture, with 
numerous avatars, claims that any G-function comes from geometry i.e. roughly 
speaking, it has an integral representation^ 

To test the validity of Christol's conjecture we look for counter-examples not 
contradicting classical conjectures. For instance, we search D-finite power series with 
integer coefficients which are not algebraic but have an integral representation and are 
not of maximal weight for the corresponding Picard-Fuchs linear ODE. 

As a first step let us limit ourselves to hypcrgeometric functions n+ iF n . The 
monodromy weight W is exactly the number of 1 among the hi. 

% The series expansion of the susceptibility of the isotropic 2-D Ising model can be recast into a series 
with integer coefficients (see [12, 16, 24, 83]), but it cannot be the diagonal of rational functions since 
the full susceptibility is not a D-finite function [83] . 
§ Bombieri-Dwork conjecture see for instance [27]. 
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3F2 



1 

-, 1 
3 



i 6 x 



1 +60 a; + 20475 a; 2 + 9373650 a; 3 



When „_i_i-F„ is globally bounded and has no integer parameters bi iW = 0), its 
minimal ODE has ap-curvature zero for almost all primes p. However, a Grothcndicck 
conjecture, proved for 3 F 2 in [85], and generalised to n +\F n in [86], asserts that, 
under these circumstances, the hypergeometric function is algebraic. We display in 
Appendix B a set of n +iF n hypergeometric functions which yield, naturally, series 
with integer coefficients, many of them corresponding to such algebraic hypergeometric 
functions. Even if such examples are quite non-trivial, the purpose of our paper is to 
focus on transcendental (non algebraic) functions. 

So we are looking for globally bounded hypergeometric functions satisfying 
1 < W < n — 1. In general such hypergeometric functions are G-series but are very 
far from being globally bounded. The hypergeometric world extends largely outside 
the world of diagonal of rational functions. 

Such an example in the first case n = 2, W = 1 was given in [43]: 

"1 4 5" 

9' 9' 9 

+ 4881796920 a; 4 + 2734407111744 a; 5 + 1605040007778900 x 6 + ••• (43) 

The integer coefficients read with the rising factorial (or Pochhammcr) symbol 

(l/9)n ■ (4/9) n ■ (5/9) n 6n p(n) 

(l/3) n .(l) B -n! p{0y [ ' 

where: 

, _ r(l/9 + n)r(4/9 + n)r(5/9 + n) 6n 

p{n) r(i/3 + n)r(i + n)r(i + n) ' 6 ■ (45j 

Note that, at first sight, it is far from clear§ on (45), or on the simple recursion on 
the p(n) coefficients (with the initial value p(0) = 1) 

P(n + l) (1 + gn) (4 + 9n) (5 + 9n) 

p{n) (l+3n)(l + n) 2 ' { ' 

to see that the p(n)'s are actually integers. A sketch of the (quite arithmetic) proof 
that the p(n)'s are actually integers, is given in [17]. 

Because of the 1/3 in the right (lower) parameters of (43), the hypergeometric 
function (43) is not an obvious Hadamard product of algebraic functions (and thus 
a diagonal of a rational function), and one can see that it is not an algebraic 
hypergeometric function either by calculating its p-curvature and finding that it is 
not zero [75], or using [87]. Proving that an algebraic function is the diagonal of a 
rational function and proving that a solution of maximal weight for a Picard-Fuchs 
equation is the diagonal of a rational function use two entirely distinct ways. The 
hope is to combine both techniques to conclude in the intermediate situation. 

This example remained for twenty years, the only "blind spot" on Christol's 
conjecture. We have recently found many other 3F2 examples:]:, such that their 
series expansions have integer coefficients but are not obviously diagonals of rational 
functions. These new hypergeometric examples are displayed in [17]. Unfortunately 
these hypergeometric examples are on the same "frustrating footing" as Christol's 
example (43): we are not able to show that one of them is actually a diagonal of a 
rational function, or, conversely, to show that one of them cannot be the diagonal of 
a rational function. 

§ In contrast with cases where binomial (and thus integers) expressions take place. 

I 2-F1 cases are straightforward, and cannot provide counterexamples to Christol's conjecture. 
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6. Integrality versus modularity: learning by examples 

A large number of examples of integrality of series-solutions comes from modular 
forms. Let us just display two such modular forms associated with HcunG functions 
of the form HeunG(a, q, 1, 1, 1, 1; x). Many more similar examples can be found in [17]. 



6.1. First modular form example 

One can, for instance, rewrite the example (37) of subsection (4.2), namely 

HeunG(-l/8, 1/4,1,1,1,1; hypergeometric function with two rational 
pullbacks: 

FeunG(-l/8,l/4,l,l,l,l;-x) = E E ( J x " ( 47 ) 

n=0 fc=0 ^ ' 

\ -1/4 

2 i RA Ji\ 



= ((1 + Ax)- (1 + 228a; + 48a; 2 + 64a; 3 ) 

/ 1 5 1728.(1 -8z) 6 -(l 

X 2^l [T77' 7o"J> 



12' 12 J! L J ' (1 + 228.x +48x 2 + 64a; 3 ) 3 • (1 +4x) 3 

1/4 



((1 -2x)-(l - 6 x + 228 x 2 -8 a; 3 )) 



, 1 5 . M1 1728- (1 - 8a;) 3 ■ (1 + xf ■ x 2 



12' 12 J! L J ' (1 -6a- + 228a; 2 -8a; 3 ) 3 ■ (1 -2a;) 3 , 

The relation between the two pullbacks, that are related by the "Atkin" involution^ 
x O — 1/8/ x, gives the modular curve: 

1953125 y 3 z 3 - 187500 y 2 z 2 ■ {y + z) + 375 y z ■ (16 z 2 - 4027 y z + 16 y 2 ) 

-64- (z + y)-(y 2 + z 2 + U87 y z) + 110592- yz = 0. (48) 

Series (47) is solution of the (exactly) self adjoint linear differential operator fl 
where (9 = x • D x ): 

x-Sl = 6 2 - x- (70 2 + 7 6 + 2) - 8a; 2 • (0 + l) 2 . (49) 



6.2. Second modular form example 

The integrality of series-solutions can be quite non-trivial like the solution of the 
Apery-like operator 

Q. = x- (1 -llx -x 2 )-Dl + {l -22x -Zx 2 )- D x - (x + 3), (50) 

or: x-n = e 2 -x- {ne 2 + ne + 3) -x 2 -(e+i) 2 , 

which can be written as a HeunG function. This (at first sight involved) HeunG 
function reads: 

TT _/ 123 55 rl/2 33 15 rl/2 i 1 1 , /ll-5 3 / 2 \ \ 

1 In previous papers [88, 28], with some abuse of language, we called such an involution an Atkin- 
Lehner involution. In fact this terminology is commonly used in the mathematical community for 
an involution r —> —N/t, on r, the ratio of periods, and not for our ^-involution. This is why we 
switch to the wording "Atkin" involution. 
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n=Ofc=0^' ^ ' ' 

but actually corresponds to a modular form, which can be written in two different 
ways using two pullbacks: 

(a; 4 + 12 a; 3 + 14 a; 2 - 12 x + 1)~ 1/4 

1728 -a; 5 ■ (1 -11; 



/,1 5. .. 



. x — x 2 ) 



(a- 4 + 12a; 3 + 14 x 2 - 12a; + l) 3 , 
(1 + 228 a; + 494 a; 2 - 228 a; 3 + a; 4 )~ 1/4 (51) 



? 4T2> 



2\5 



5 , r , 1728 -a: - (1 - 11 x - x 2 ) 

x 2-F1 ' 1 



(1 + 228 a- + 494 x 2 - 228 x 3 + a; 4 ) 3 



More modular form examples of series with integer coefficients are given in [17]. 
The modular form examples displayed in Appendix A of [17], correspond to lattice 
Green functions [46]. Therefore, they have n-fold integral representations, and, after 
section (3.4), can be seen to be diagonals of rational functions. In contrast the 
modular form examples displayed in Appendix H of [17] correspond to differential 
geometry examples discovered by Golyshev and Sticnstra [89] , where no n-fold integral 
representation is available at first sight. 



7. Integrality versus modularity 

7.1. Diffeomorphisms of unity pullbacks 

Let us consider a first simple example of a hypergeometric function which is 
solution of a Calabi-Yau ODE, and which occurred, at least two times in the 
study of the Ising susceptibility n-fold integrals [28, 29] x^ an d Xd > namer y 
4F3QI/2, 1/2, 1/2, 1/2], [1, 1, l];256x), where we perform a (diffcomorphism of 
unity) pullback: 

*F 3 ([llLk,[hl,l}; t- 2 f X 2 - ) = 1 +16-* (52) 

V 2 2 2 2 1 + c\ x + C2 x z + ■ ■ ■ / 

+ 16 -(81- c^-a; 2 + 16 • (10000 + c\ - c 2 - 162 Ci ) • a; 3 +••• 

If the pull-back in (52) is such that the coefficients c„, at its denominator, are integers, 
one finds that the series expansion is actually a series with integer coefficients, for every 
such pullback (i.e. for every integer coefficients c„). Furthermore, a straightforward 
calculation of the corresponding name q(x) and its compositional inverse (mirror map) 
x(q), also yields series with integer coefficients: 

q (x) = x +(64-ci)-a; 2 +(c 2 + 7072 -c 2 -128ci)-a; 3 +■••, (53) 

x(q) = q + (ci - 64) • q 2 + (c 2 + 1120 + c 2 - 128 a) ■ q 3 + ■■■ , (54) 

when its Yukawa coupling [28], seen as a function of the nome q, K(q) is also a series 
with integer coefficients and is independent of the pullback: 

K(q) = 1 + 32 ■ q + 4896 • q 2 + 702464 • q 3 + ■■■ (55) 

This independence of the Yukawa coupling with regards to pullbacks, is a known 
property, and has been proven in [70] , for any pullback of the diffcomorphism of unity 
form p(x) = x + • • ■ 
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Seeking for Calabi-Yau ODEs, Almkvist et al. have obtained [73] a quite large list 
of fourth order ODEs, which are MUM by definition and have, by construction, the 
integrality for the solution-scries analytic at x = 0. Looking at the Yukawa coupling 
of these ODEs is a way to define equivalence classes up to pullbacks of ODEs sharing 
the same Yukawa coupling. This "wraps in the same bag" all the linear ODEs that are 
the same up to pullbacks. Let us recall how difficult it is to see if a given Calabi-Yau 
ODE has, up to operator equivalence, and up to pullback, a hypcrgcomctric function 
solution [28, 29], because finding the pullback is extremely difficult [28, 29]. We may 
have, for the Ising model, some n +iF n hypergeometric function prejudice [28, 29]: it 
is, then, important to have an invariant that is independent of this pullback that we 
cannot find most of the time. 

Finally, let us remark that the Yukawa coupling is not preserved by the operator 
equivalence. Two linear differential operators, that arc homomorphic, do not 
necessarily have the same Yukawa coupling (see Appendix C). 



7.2. Yukawa couplings in terms of determinants 

Another way to understand this fundamental pullback invariance, amounts to rewriting 
the Yukawa coupling [70, 90], not from the definition usually given in the literature 
(second derivative with respect to the ratio of periods), but in terms of determinants 
of solutions (Wronskians, ...) that naturally present nice covariance properties with 
respect to pullback transformations (see Appendix C). 

We have the alternative definition for the Yukawa coupling given in Appendix C: 

*<•> - - ^ 

where the determinantal variables W m 's are determinants built from the four solutions 
of the MUM differential operator. This alternative definition, in terms of these 
Wm's, enables to understand the remarkable invariance of the Yukawa coupling by 
pullback transformations [29] . These determinantal variables W m quite naturally, and 
canonically, yield to introduce another "Yukawa coupling" (which, in fact, corresponds 
to the Yukawa coupling of the adjoint operator (see C.12)). This "adjoint Yukawa 
coupling" is also invariant by pullbacks. It has, for the previous example, the following 
scries expansion with integer coefficients: 

K*(q) = I + 32 ■ q + 4896 ■ q 2 + 702464 • q 3 +■■■ (57) 

which actually identifies with (55). The equality of the Yukawa coupling for this order- 
four operator, and for its (formal) adjoint operator, is a straightforward consequence 
of the fact that the order-four operator annihilating 4-F 3 ([i, |, |, j], [I, 1, I]; 256x) 
is exactly self-adjoint, and, more generally, of the fact that the order-four operator, 
annihilating (52), is conjugated to its adjoint by a simple function. 



7.3. Modularity 

This example, with its corresponding relations (52), (53), (55), (57) may suggest 
a quite wrong prejudice that the integrality of the solution of an order-four linear 
differential operator automatically yields to the integrality of the nome, mirror map 
and Yukawa coupling, that we will call, for short, "modularity. This is far from being 
the case, as can be seen, for instance, in the following interesting example, where the 
nome and Yukawa coupling K(q) do not correspond to globally bounded series, when 
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the 4^3 solution of the order-four operator as well as the Yukawa coupling seen as a 
function of x, K(x), are, actually, both series with integer coefficients. 

Let us consider the following 4^3 hypergeometric function which is clearly a 
Hadamard product of algebraic functions and, thus, the diagonal of a rational function: 

Mih Z> I'M,!];*) 

= (1 - xr 1 / 3 * (1 - x)- 1 ' 2 * (1 - x)- 1 ^ * (1 - x)- 3 / 4 

= Diag((l - Zl )-^ (1 - z 2 )-Va (1 _ ^-V* (1 _ z 4 )~ 3 / 4 ) , 
It is therefore globally bounded: 
4 F 3 ([i, i, -, [1, 1, 1]; 2304a;) = 1 +72a; + 45360a; 2 + 46569600a; 3 

+ 59594535000 a; 4 +86482063571904 a; 5 +••• (58) 

Its Yukawa coupling, seen as a function of x, is actually a series with integer 
coefficients in a;: 

= 1 +480 a; + 872496 x 2 + 1728211968 a; 3 + 3566216754432 a; 4 

+ 7536580798814208 a; 5 + 16177041308360579328 a; 6 + ••• (59) 

However, do note that the series, in term of the nome, is not globally bounded: 

K(q) = 1 + 480 q + 653616 q 2 + 942915456 q 3 + 1408019875200 q 4 + 

12 

+ 571436303929319146711343817202689132288 -|— + • • • (60) 

In fact, the nome q(x), and the mirror map x(q), are also not globally bounded. 
Note that in this example, the non integrality appears at order twelve (for x(q), q(x) 
and K(q)). If the prime 11 in the denominator in (60) was the only one, one could 
recast the series into a series with integer coefficients introducing another rescaling 
2304 a; — > 11 X 2304 x. But, in fact, we do see the appearance of an infinite number 
of other primes at higher orders denominators in x(q), q(x) and K{q). 



7.4- Order-two differential operators uj n associated with modular forms 

After Maier [91] let us underline that modular forms can be written as hypergeometric 
functions with two different pullbacks, and, consequently, one can associate order-two 
differential operators to these modular forms. 
Let us consider the two order-two operators 

„ 2 . Dl + ^LtlL. Dl +75 4 Tr , (6i) 

(64 x + 1) • x (64 x + 1) x 

-3= Dl + r9 ( 7 45 ^+ 1) D x + 3 , (62) 
(27x + l)-x (27a; + l)a; 

which are associated with two modular forms corresponding, on their associated nomcs 
q, to the transformations q — > q 2 and q — > q 3 respectively (multiplication of r, the 
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ratio of their periods by 2 and 3), as can be seen on their respective solutions: 
r l 1, r „, i , _ 1M „ / r 1 5, N1 1728a; 



[l];-64x) = (l+256x)- 1 / 4 - 2 F 1 ([^, A], [1]; 



and: 



4' 4 J ' L J ' ' K ' V12' 12 J! L J ' (1 + 256x) 

= ( 1 + w ')- 1/4 -^(i5-s"-W'(r^) (63) 

1 -4a; +100 a; 2 - 3600a; 3 + 152100 x 4 - 7033104a; 5 + 344622096 x 6 
- 17582760000 x 7 + 924193822500 a; 8 - 49701090010000 a; 9 + ■ ■ ■ 



Via „ /, 1 5, N1 1728a; 



((1 + 27x)(l + 243x) 3 ) • 2 F 1 ([^, A], 



(1 + 243x) 3 (l + 27x) 



\ -1/12 / 1 ^ 1 798 t 3 \ 

(l + 27*)(l + 3*)») ■ 2Fl ([i JL], [1]; Ii - F ^£_) (64) 



= 2^1 ( ^], [1], -27 x) = 1 -3a; +36a^ - 588a: 3 + 11025 a; 4 - 223587a; 5 

+ 4769856 a; 6 - 105423552 a; 7 + 2391796836 a; 8 - 55365667500 a; 9 
The relation between the two Hauptmodul pullbacks in (63) 
1728 x 1728 x 2 / 1 



(65) 



(1 + 256x) 3 ' (1 +16x) 3 

corresponds to the (genus-zero) fundamental modular curve: 

5 9 • uV - 12 • 5 6 • u 2 v 2 ■ (u + v) + 375 uv ■ (16 u 2 + 16 v 2 - 4027 uv) 

-6A(u+v)- (v 2 + 1487 uv + u 2 ) + 2 12 3 3 -uv = 0. (66) 
The relation between the two Hauptmodul pullbacks in (64) 

1728 x 1728 x 3 / 1 \ 

U ~ (1 +243x) 3 (l +27 x)' V (1 +3 a;) 3 (l +27 x) U V3 6 x/' ( ' 
corresponds to the (genus-zero) modular curve: 
2 27 5 9 • u 3 v 3 ■ (u+v) + 2 18 5 6 u 2 v 2 ■ (27 v 2 + 27 u 2 - 45946 to) 

+ 2 9 3 5 5 3 uv (u+v) - (v 2 + 241433 uv + u 2 ) (68) 
+ 729 (u 4 + u 4 ) - 39628 • 3 9 • (u 2 + v 2 ) ■ uv + 15974803 ■ 2 • 3 10 • u 2 v 2 
+ 31 • 2 9 3 n uv ■ (u +v) -2 12 3 12 uv = 0. 

Similarly, one can consider the order-two operators w n associated with other 
modular forms corresponding to r — s- n ■ r. The w n 's can be simply deduced from 
Maier [91], for modular forms corresponding to genus-zero curves i.e. for n = 2, 3, 
4, 5, 6, 7, 8, 9, 10, 12, 13, 16, 18, 25. Since the solutions can be written as 2-F1 
hypcrgcometric up to rational pullbacks, these genus-zero w n 's are obviously order- 
two operators. After a simple rescaling, the solutions analytic at x = 0, can be 
rewritten as series with integer coefficients. 

One can also consider the other w„'s corresponding to higher-genus modular 
curves. In these cases, one does not have a rational parametrisation like (65) or (67), 
but one still has an identity of the same hypergeometric function with two different 
pullbacks, these two pullbacks being algebraic functions and not rational functions 
(see (65) or (67)). These algebraic functions correspond to the so-called modular 
polynomials [17]. 
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For instance for r — > 11 • r, one has a genus-one modular curve, the modular 
polynomial reads: 

P^ (x,H) = (1 + 228 x + 486 x 2 - 540 x 3 + 225 a; 4 ) 3 • H 2 (69) 
- 1728 • Qx(x) ■ x ■ H + 1728V 2 , with: 

Qi{x) = 1 -55a; + 1188a; 2 -12716a; 3 + 69630a; 4 - 177408 x 5 + 133056 x 6 
+ 132066 a; 7 - 187407 a; 8 +40095 a; 9 + 24300 a; 10 -6750X 11 . 

One has the identity 

2Fl ( [ iT il ] ' [1];Hl ) = ^•^([^.^l.W;^) 5 (70) 

where the two Hauptmoduls i?i and #2 are the two solutions of P*i(x, H) = 0, 
and where A(x) is the algebraic function such that 

A(x) 4 ll 2 2 7321 - 87612 a; + 73206 a; 2 + 21060 a; 3 - 23175 a; 4 

ll 2 + A(x) 4 Tl 1 ' 1 + 228 a; + 486 x 2 - 540 a; 3 + 225 x 4 ' 

this last relation between A(x) 4 and x corresponding to a genus-one curve with the 
same j -invariant as the genus-one curve Pii(x, H) = 0, namely j = — 496 3 /ll 5 . 
The two hypergeometric functions in (70) are actually series with integer coefficients: 

2 Fi([l/12, 5/12], [1]: #1) = 1 +60.T -4560a; 2 + 614400a; 3 - 95660400x 4 

+ 16231863060a; 5 - 2905028387700a; 6 + ••• (71) 

2 .Fi([l/12, 5/12], [1]; H 2 ) = 1 + 60 a; 11 + 3300 a; 12 + 110220a; 13 + 2904660 x 14 
+ 66599940 a; 15 + 1394683620 a; 16 + 27425371380 x 17 + ••■ (72) 

More details are given for this t — > 11 • r case in Appendix I of [17]. Note that 
although 2-Fi([l/12, 5/12], [1]; H x ) and 2^1 ([1/12, 5/12], [1]; H 2 ) are solutions of the 
same order- four operator [17], one can find an appropriate algebraic function A(x), 
such that .4(a;) • 2-Pi([l/12, 5/12], [1], Hi) is solution of an order-two operator u)u 
(see [17] for more details). 

The other w„'s, corresponding to higher genus modular curves [92], are actually 
also order-two operators. The explicit expressions of w„'s for the elliptic values 
n = 17, 19, and the hyperelliptic values [92] n — 23, 31, 39, 41, 47, 59, 71 are given 
in [17]. The genus of the associated modular curves [92], is respectively [17] genus- 
one for u>n(x), cDig(x), genus-two for £23(3;), 0)29(3;), 0)31(2;), genus-three for 0)41(2;), 
genus-four for 0)47(2;), genus-five for 0)59 (x), and genus-six for 0)71(3;). 

7. 5. Hadamard products of u n 's 

The two operators 0J2 and 0J3 have a "modularity" property: their series expansions 
analytic at x = 0, (63) and (64), as well as the corresponding nomcs, mirror maps are 
series with integer coefficients. The Hadamard product is a quite natural operation to 
introduce because it preserves the global nilpotence of the operators, it preserves the 
integrality of series-solutions, and it is a natural operation to introduce when seeking 
for diagonals of rational functions^. Let us perform the Hadamard product of these 



And, consequently, has been heavily used to build Calabi-Yau-like ODEs (see Almkvist et al. [70]). 
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two operators. With some abuse of language [29], the Hadamard product of the two 
order-two operators (61) and (62) 

H { j + t) n 4 | a (2064£-l) 3 (19020 x- 7) 2 

(4788 a; -1) 12 . , 

D * + TT^Z 71 ^> ( 73 ) 



(1728 a;- 1) - a; 3 (1728 a;- 1) 

is defined as the (minimal order) linear differential operator having, as a solution, the 
Hadamard product of the solution-series (63) and (64), which is, by construction, a 
series with integer coefficients. This series is, of course, nothing but the expansion of 
the hypcrgcomctric function: 

|, |, |], [1,1,1]; 1728 x) (74) 



2 F 1 ([\, J], [1]; -64 re) * 2 ^([i 1], [1]; -27a;). 



In a similar way one can consider (see [17]) H2.2 = ^2 * ^2 (resp. .2/3,3 = 
W3 7k- W3) the Hadamard product of the order-two operator (61) (resp. (62)) with 
itself (Hadamard square). These two operators have respectively the hypcrgcometric 
solutions 

* F ^\> i> i> t 1 - ^ 4096 *)> ^([^ s> 5> 5]' t 1 - x > i]; 729a; ). (75) 

corresponding to scries expansions with integer coefficients. These operators H2.2, 
H3 3 arc MUM operators. We can, therefore, define, without any ambiguity, the nomc 
(and mirror map) and Yukawa coupling of this order- four operator [29] . One finds out 
that the nomef, and the mirror map (and the Yukawa coupling as a function of the x 
variable), are not globally bounded: they cannothe reduced, by one rescaling, to series 
with integer coefficients. 

The three linear differential operators i/2,3, i/2,2 and , are MUM and of 

order four, however, they are not of the Calabi-Yau type. 

7. 6. Hadamard products versus Calabi- Yau ODEs 

The occurrence of Calabi-Yau type operators, that we could imagine, at first sight, 
to be extremely rare, is in fact quite frequent among such Hadamard products, 
as can be seen with other values of n and m. For instance, one can introduce^ 
i?4.4 = -k LU4. , the Hadamard square of W4, which is an irreducible order-four 
linear differential operator, and has the hypergeometric solution already encountered 
for some n-fold integrals of the decomposition of the full magnetic susceptibility of 
the Ising model [28, 29] (see also subsections (7.1) and (7.2)): 

1. I];256a0 (76) 



L 2' 2' 2' 2 J 



[1]; -16*) * \], [1]; -16s). 



f The nome of the Hadamard product of two operators has no simple relation with the nomc of these 
two linear differential operators. 

\ To get the Hadamard product of two linear differential operators use, for instance, Maple's 
command gfun[hadamardproduct]. 
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The associated operator having (76) as a solution, obeys the "Calabi-Yau 
condition" that its exterior square is of order five. 

Let us give in a table the orders (which go from 4 to 20) of the various 
#m,n = #n,m Hadamard products of the order-two operators associated with the 
(genus-zero) modular forms operators lj„ and uj m : 
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where the star * denotes Calabi-Yau ODEsjj. 

The following operators arc of order four: i?2,2, #2,3; #2.4, #2, 6: #2,8> #2,9) 
#3,31 #3,4 1 #3,6, #3.8, #3,9; ■•■ Their exterior squares, which are of order six, do not 
have rational solutionsf . 

The order-four operators #3,3, #3,4, are all MUM operators^, but, similarly to 
the situation encountered with ^2,2, their nome, mirror map and Yukawa couplings 
are not globally bounded. 

The following operators are of order six: i?2,5, #2,7, #3.5, #3,7, #4,5, #4,7, 
#5,5; #5,6> #5,8: #5.9, #6,7: #7,7> #7,8: #7,9: •■• Their exterior square, which are 
of order fifteen, do not have rational solutions (and cannot be homomorphic to higher 
order Calabi-Yau linear ODEs). 

Remarkably the following ten order- four operators i?4.4, #4,6: #4,8; #4.9) #6.6 , 
#6,8, #6,9, #8,8, #8,9, #9,9 (with a star in the previous table) are all MUM, and are 
such that their exterior squares are of order five^: they are Calabi-Yau ODEs. Actually 
the nome, mirror map and Yukawa coupling series are series with integer coefficients 
for all these order-four Calabi- Yau operators. Their Yukawa coupling and their adjoint 
Yukawa coupling identify. The Yukawa coupling series of these Calabi-Yau operators 
are respectively, for #4,4 

K(q) = K*(q) = 1 + 32 • q +4896- q 2 + 702464 • q 3 + 102820640 ■ g 4 

+ 15296748032 • q 5 + 2302235670528 • q 6 + ■■■ (77) 

t) Recall that Calabi-Yau ODEs are defined by a list of constraints [70], the most important ones 
being, besides being MUM, that their exterior square are of order five. There are more exotic 
conditions like the cyclotomic condition on the monodromy at 00, see Proposition 3 in [73]. 
f They cannot be homomorphic to Calabi-Yau ODEs. 

Note that the Hadamard product of two MUM ODEs is not necessarily a MUM ODE: the order-six 
operator H3 7 is not MUM. 

§ They are conjugated to their (formal) adjoint by a function. 
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which is #3 in Almkvist et al. large tables of Calabi-Yau ODEs [73], and is the 
well-known one for 4-^3 ([1/2, 1/2, 1/2, 1/2], [1, 1, 1]; 256 a;), and for F 4 ,6 

K(q) = K*{q) = 1 + 20 -q+ 36 • q 2 + 15176 • q 3 + 486564 • q 4 

+ 21684020 • q 5 + 1209684456 • q 6 + ■■■ (78) 

which is #137 in tables [73]. 

We give, in [17], the expansion of the Yukawa coupling for a set of other H m , n 
such that their exterior squares are order five (not six as one could expect for a generic 
irreducible order- four operator), that actually are Calabi-Yau operators. Actually 
operator H^s is #36 in Almkvist et al. large tables of Calabi-Yau ODEs [73]. 
Operators i/4,9, H 6t6 , i^g and H e _ 9 are respectively [73] #133, #144, #176 and 
#178. Furthermore, operators -£/g,g, -Hg,g and Hg^ are respectively [73] #107, #163 
and #165. 

It will be shown, in a forthcoming publication, that the occurrence of an order five 
for the exterior power (the "Calabi-Yau condition") means that these operators are 
necessarily conjugated (by an algebraic function) to their adjoints. Thus, the "adjoint 
Yukawa coupling" K*{q) is necessarily equal to the Yukawa coupling K(q) for these 
operators. 

On the other hand, the ten linear differential operators denoted by a star * in the 
previous table all share the same property: they have, as a solution, the Hadamard 
product of two HeunG functions solutions of the form HeunG(a, q, 1, 1, 1, 1; x). 
Note, however, that this HeunG-viewpoint of the most interesting if mi „'s does not 
really help. Even inside this restricted set of HeunG functions solutions of the form 
HeunG(a, q, 1, 1, 1, 1; x) it is hard to find exhaustively the values of the parameter a 
and of the accessory parameter q, such the series HeunG(a, q, 1, 1, 1, 1; x) is globally 
bounded, or, just, such that the order- two operator, having HeunG(a, q, 1, 1, 1, 1; x) 
as a solution, is globally nilpotent [25]. 

Many H m n are not MUM, for instance the order-eight operator -£/i2,i2, or the 
order-six operator ££3,7, are not MUM. Concerning -^3,7, and as far as its six solutions 
are concerned, they are structured "like" the four solutions of an order-four MUM 
operator, together with the two solutions of another order-two MUM operator, but the 
order-six operator -££3.7 is not a direct-sum of an order-four and order-two operator. 
We have two solutions analytical at x = (with no logarithmic terms), and two 
solutions involving ln(cc). A linear combination of these two solutions analytical at 
x = is, by construction a series with integer coefficients (the Hadamard product 
of the two series with integer coefficients which are the initial ingredients in this 
calculation), when the other linear combinations are not globally bounded. 

8. Calabi-Yau Modularity 

The previous examples correspond to a "modularity" inherited from elliptic curves, 
more precisely Hadamard products of modular forms. Let us consider, here, two 
Calabi-Yau examples that do not seem to be reducible! to 4F3 hypergeometric 
functions. 

\ The possibility that a solution of an order-four operator, non-trivially equivalent to these Calabi- 
Yau operators [72], could be written as a 4.F3 hypergeometric function, up to an involved algebraic 
pullback, is not totally excluded. However it is extremely difficult to rule out such highly non-trivial 
hypergeometric scenario. 
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A first order- four Calabi-Yau operator, found by Batyrev and van Straten [72], 
which is self-adjoint and also corresponds to Hadamard products of simple 
hypergeometric functions (see (D.2)), is given in Appendix D. All the associated series 
(solution (D.2), nome, Yukawa coupling) are series expansion with integer coefficients. 
We do not have a representation of the solution (D.2), as an n-fold integral of the 
form (34). However, since (D.2) can be expressed as a sum of products of binomials 
(see (D.3)), we can conclude, again, that (D.2) is actually a diagonal of a rational 
function. 

8.1. A Batyrev and van Straten Calabi-Yau ODE [72] 

A second example of order-four operator, corresponding to Calabi-Yau 3-folds in 
Pi x Fi x ?! x has been found by Batyrev and van Straten [72] (seef page 34): 

B 2 = 4 -Ax- (5 2 + 50 +2) • (2 + l) 2 

+ 64a; 2 ■ (20 + 3) ■ (20 + 1) ■ (20 + 2) 2 . (79) 

It corresponds to the series-solution with coefficients: 



(2n\ A fn\ 2 (2k\ (2n-2k 



x X fe = xx XXX X 

Its Wronskian Wa is a rational function such that: 

* 3 ' = (1-64,) 1 (1-16,) - (81) 

This operator is also a Calabi-Yau operator: it is MUM, and it is such that its exterior 
square is order five. This order five property is a consequence of B 2 being conjugated 
to its adjoint: B 2 ■ x — x ■ adjoint(B 2 ). 

The series-solution of (79) can be written as an Hadamard product 

S = (1 -4.t)~ 1/2 * HeunG(4, 1/2,1/2,1/2,1,1/2; 16a;) 2 (82) 
= 1 + 8 x + 168 x 2 + 5120 x 3 + 190120 x 4 + 7939008 x 5 + 357713664 x 6 + • • • , 

the modular form character of HeunG(A, 1/2, 1/2, 1/2, 1, 1/2; 16 x) being illustrated 
with identities (A. 4) in Appendix A of [17]. Its nome reads: 

q = x + 20 x 2 + 578 x 3 + 20504 x 4 + 826239 x 5 + 36224028 x 6 + 1684499774 x 7 

+ 81788693064 a; 8 + 4104050140803 a; 9 + 211343780948764 x 10 + ■■■ (83) 

The mirror map of (79) reads: 

x(q) = q -20q 2 + 222 q 3 -21Mq 4 + 21293 q 5 - 307224 q 6 + 80402 q 7 

- 67101504 q s - 1187407098 q 9 - 37993761412 q w + ■■■ (84) 

The Yukawa coupling of (79) reads: 

K(q) = K*(q) = 1 +4g +164q 2 + 5800 q 3 + 196772 q 4 + 6564004 q 5 

+ 222025448 q 6 + 7574684408 q 7 + 259866960036 q 8 + ••• (85) 

The equality of the Yukawa coupling with the "adjoint" Yukawa coupling, K(q) = 
K*(q), is a straight consequence of relation B 2 ■ x = x ■ adjoint(B 2 ). 

f There is a small misprint in [72] page 34: (2 8 + 1) must be replaced by (2 9 + l) 2 in the Ax term. 
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Recalling Batyrev and van Stratcn [72], and following Morrison [71], do note that 
one can also write the Yukawa coupling as: 

= xjgf -Wl'\ { g dx^ = <!.(,. ^l) 3 , (86) 
Vo V ^(?) d q J Vo V dq J 

where Wa is the Wronskian (81). From this alternative expression for the Yukawa 
coupling, valid when the operator is conjugated to its adjoint (see (C.7)), it is obvious 
that if the analytic series ya(x), as well as the nome (83) are series with integer 
coefficients, then, the mirror map (84) is also a series with integer coefficients, and, 
therefore, yo seen as a function of the nome q, as well as x 3 W^ 2 (since it is a rational 
function) are also series with integer coefficients. Consequently, the Yukawa coupling 
is a series with integer coefficients (as a series in q or in x). 

More generally, if one assumes that a linear differential operator has a globally 
bounded solution-series, one knows that this operator is a G-operator, necessarily 
globally nilpotent, and, consequently, its Wronskian, or the square root of the 

1/2 

Wronskian (see W 4 in (86)), will be a A^-th root of a rational function, and, thus, 
will correspond to a globally bounded series. Thus, the globally bounded character of 
the analytic series yo(x) together with the nome, yields the globally bounded character 
of the mirror map, Yukawa coupling, that we associate with the modularity. Similar 
results can be found in Delaygue's thesis [93], in a framework where the coefficients 
of hypcrgcomctric series are ratio of factorials (see (Appendix B)). 

In contrast the globally bounded character of the analytic series yo(x), together 
with the globally bounded character of the Yukawa coupling (seen for instance as a 
series in x) does not imply that the nome, or the mirror map, are globally bounded 
as can be seen on example (58) (see (59) and (60)). 



8.2. An operator non trivially homomorphic to B 2 
Let us, now, consider the order-four operator 
B 2 = 256a; 2 - 9 2 (2 9 + 3) (2 9 + 1) 

-Ax- (26> + l) (26»-l) (56» 2 - 59 + 2) + (9 - l) 4 . (87) 
This operator is non-trivially| homomorphic to the Calabi-Yau operator (79): 

B 2 -x-(20+l) = x-(29 +!)■ B 2 . (88) 

As a consequence of the previous intertwining relation, one immediately finds that the 
series-solution analytic at x = of this new MUM operator (87) is nothing but the 
action of the order-one operator x ■ (29 + 1) on the series (82), and reads: 

x- (29 = x + 24a; 2 +840 x 3 +35840 x 4 + 1711080a; 5 +87329088 x 6 

+ 4650277632 x 7 + 254905896960 a; 8 + ••• (89) 

It is obviously also a scries with integer coefficients (the action of x ■ (2 9 + 1) on the 
series with integer coefficients is straightforwardly a series with integer coefficients). 
More generally, the globally bounded series remain globally bounded series by non 
trivial operator equivalence, namely homomorphisms between operators (generically 
the intcrtwincr operators are not simple functions). 

The exterior square of the order-four operator (87) is an order-six operator which 
is, in fact, the direct sum of an order-five operator £5 and an order-one operator. 



I The intertwiners between B2 and B2 are operators not simple functions. 
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Operator B 2 is non-trivially homomorphic to its adjoint: 

B 2 ■ x 3 ■ (2 9 + 3) • (2 6» +5) = x 3 ■ (2 6 + 3) • (2 6 + 5) • adjoint(B 2 ). (90) 

The Yukawa coupling of this order-four operator (87), non-trivially homomorphic to 
(79), reads: 

a 4 a 5 

K{q) = 1 -4q -UOq 2 -AOAOq 3 - 64436 -^- + 1889332 ^- 

+ 88331368 — + 1652707624 69295027684 — + • • • (91) 

5 9 63 

The Yukawa coupling scries (91) is not globally bounded. 

The "adjoint Yukawa coupling" of this order-four operator (87) reads: 

K*(q) = 1 +12 a + 564g 2 + 20440 q 3 + 865732 q 4 + 37162444 q 5 

+ 8255346664 — + 1121762648248 %- + 72336859374772 %- + ■■■ (92) 
5 15 21 

Again, the adjoint Yukawa coupling series (92) is not globally bounded. 

On this example one sees that the Yukawa coupling of two non-trivially 

homomorphic operators are not necessarily equal. The Yukawa couplings of two 

homomorphic operators are equal when the two operators are conjugated by a function 

(trivial homomorphism) . The modularity property is not preserved by (non-trivial) 

operator equivalence: it may depend on a condition that the exterior square of the 

order-four operators is of order five. The Calabi-Yau property is not preserved by 

operator equivalence. 

To sum-up: All these examples show that the integrality (globally bounded 
scries) is far from identifying with modularity. 

9. Conclusion 

Seeking for the linear differential operators for the X s i we discovered, some 
years ago, that they were Fuchsian operators [11, 13], and, in fact, "special" 
Fuchsian operators, namely Fuchsian operators with rational exponents for all their 
singularities, and with Wronskians that are N-th roots of rational functions. Then 
we discovered that they were G-operators (or equivalently globally nilpotent [25]), 
and more recently, we accumulated results [29] indicating that they are "special" G- 
operators. There are, in fact, two quite different kinds of "special features" of these 
G-operators. On one side, we have the fact that one of their solutions is not only 
G-series, but is a globally bounded scries. This special character has been addressed 
in this very paper, and we have seen that, in fact, this "integrality property [94]" is a 
consequence of quite general mathematical assumptions often satisfied in physics (the 
integrand is not only algebraic but analytic in all the variables (35)). However, we 
have also seen another special property of these G-operators, namely the fact that 
they seem to be quite systematically homomorphic to their adjoints [29] . We will show, 
in a forthcoming publication, that this last property amounts, on the associated linear 
differential systems, to having special differential Galois groups, and that their exterior 
or symmetric squares, have rational solutions. This last property is a property of a 
more "physical" nature than the previous one, related to an underlying Hamiltonian 
structure [95], or as this is the case, for instance in the Ising model, related to the 
underlying isomonodromic structure in the problem, which yields the occurrence of 
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some underlying Hamiltonian structure [95] . In general the integrality of G-operators 
does not imply the operator to be homomorphic to its adjoint, and conversely being 
homomorphic to its adjoint does not imply:}: integrality (and even does not implyf 
the operator to be Fuchsian). Interestingly, the X s : as weu as many important 
problems of theoretical physics, correspond to G-operators that present these two 
complementary "special characters" (integrality and, up to homomorphisms, self- 
adjointness), and, quite often, this is seen in the framework of the emergence of 
"modularity" . 

Nomes, mirror maps, and Yukawa couplings are not D-finite functions: they are 
solutions of quite involved non-linear (higher order Schwarzian) ODEs (see for instance 
Appendix D in [25]). Therefore, the question of the series integrality of the nomes, 
mirror maps, Yukawa couplings, and other pullback-invariants (see Appendix C) 
requires to address the very difficult question of series-integrality for (involved) non- 
linear ODEs. Note, however, as seen in Section 8.1, in particular in (86), that the 
integrality of the series yo(x) and of the nome q(x) are sufficient to ensure, provided the 
operator is conjugated to its adjoint (see (C.7)), the integrality of the other quantities 
such as the Yukawa coupling, mirror maps. However the integrality of the nome 
remains an involved problem. These questions will certainly remain open for some 
time. 

In contrast, and more modestly, we have shown that a very large sets of problems 
in mathematical physics (see sections (3.4), (4) and (5.1)) actually correspond to 
diagonals of rational functions. In particular, we have been able to show that the 
-fold integrals of the susceptibility of the two-dimensional Ising model are 
actually diagonals of rational functions for any value of the integer n, thus proving 
that the x^'s arc globally bounded for any value of the integer n. As can be seen 
in the "ingredients" of our simple demonstration (see (3.4)), no elliptic curves, and 
their modular forms [97], no Calabi-Yau [98], or Frobcnius manifolds [95], or Shimura 
curves, or arithmetic lattice assumption [99, 100] is required to prove the result. We 
just need to have an n-fold integral such that its integrand is not only algebraic, but 
analytic in all the variables. 
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Appendix A. Seeking for the "minimal" rational function 

The effective calculations of section (4.2) guarantee to obtain an explicit expression for 
the rational function associated with (40), however the rational function is far from 
being unique. Recalling the well-known Apery series A(x), and its rewriting due to 

| See Appendix M and Appendix O in [17] which give an example of a (hypcrgcomctric) family of 
order-four operators satisfying the Calabi-Yau condition that their exterior square is of order five, and, 
even, a family of self-adjoint order-four operators, the corresponding hypcrgcomctric solution-series 
being not globally bounded. 

f For instance the operator D™ — x D x — 1/2 (see page 74 of [96] ) with an irregular singularity is 
self- adjoint. 
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Strehl and Schmidt [80, 81, 82], 

°° n /r>\ 2 / n +A 2 oo n k , s / . , , 3 

^) = EE © ( ; ) = EEE (;)( r 

n=0 A;=0 V ' V 7 n=0 k=0 j=0 ^ ' V 

1 +5a; + 73 a: 2 + 1445a; 3 + 33001 a; 4 +•■•, (A.l) 

_4(a;) is known to be the diagonal of the rational function in five variables I/R1/R2 
where R\, R 2 read [55]: 

Rl = 1 - Z , R.2 = (1 - Zl)(l - Z 2 )(l - 23) (1 - Z 4 ) - ZoZxZ 2 , 

as well as the diagonal of the rational function in five variables I/Q1/Q2 where Qi, Q 2 
read [56, 43]: 

Qi = 1 - z 1 z 2 z 3 z 4 , Q 2 = (1 - z 3 )(l - z 4 ) - z - (1 + «i)(l + z 2 ), 

and aZso the diagonal of the rational function in six variables 1 / P\ / P 2 / P3 where 
P u P 2l P 3 read [55]: 

Pi = 1 — Zq Zi, P% = 1 — z 2 — z 3 — z z 2 z 3 , P 3 = 1 — z 4 — 25 — Z\ z 4 z 5 . 

A yet different diagonal representation for the Apery series, due to Delayguef, is 
provided by the diagonal of the rational function in eight variables: 

1 

(1 - z 4 z 5 z 6 z 7 ) ■ (1 - z ■ (1 + z A )) ■ (1 - z x ■ (1 + z 5 )) ■ (1 - z 2 - z 6 ) ■ (1 - z 3 - z 7 ) ' 

Calculations similar to (39) on these new binomial expressions provide two new 
rational functions such that (A.l) can be written as the diagonal of one of these two 
rational functions. One is a rational function of five variables, of the form 1 / /Q 2 

<2i 5) = 1 • (1 + Zl) (1 + Z 2 ) (1 + Z3) (1 + z 4 ), 

Q { 2 ] = 1 -zo ■ (l + zi)(l + Z2)(l + ^) 2 (l + ^) 2 , (A.2) 
and the other one, is a rational function of six variables, of the form 1/Q\ /Q 2 /Q 3 
Q { i } = 1 - z z 3 z 4 z 5 ■ (1 + zi) (1 + z 2 f (1 + 23) (1 + z 4 ) (1 + z 5 ), 

Q 2 6) = 1 - z Z 4 Z 2 Z 3 Z 4 Z 5 ■ (1 + Zi) (1 + z 2 ), 

Q { 3 } = 1 - zo ■ (1 + z x ) (1 + z 2 f (1 + z 3 ) (1 + 24) (1 + 25). (A.3) 

Appendix B. Hypergeometric series with coefficients ratio of factorials 

As a consequence of the classification by Beukers and Heckman [86] of all algebraic 
n F n —is, the 8-F7 hypergeometric scries 



^7 



1 2 3 4 1 2 1 

5' 5' 5' 5' 2' 3' 3 



2 14 3 9 5 5 x 



1 7 11 13 17 19 23 29 

30' 30' 30' 30' 30' 30' 30' 30 
has integer coefficients, and is an algebraic function. The Galois group belonging to 
this function is the Weyl group W(E$) which has 696729600 elements [101]. It is 
an algebraic series of degree 483840. More precisely, it was noticed by Rodrigucz- 
Villcgas [102] that the previous power series reads: 

V (30n)! " ! x n (B 1) 

^ (15n)!(10n)!(6n)l ' 1 ' ' 

f Personal communication. 
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which is closely related to the series introduced by Chebyshev in his work [103] on 
the distribution of prime numbers to establish the estimate [17] on the prime counting 
function tt(x). 

Considering hypergeometric series such that their coefficients are ratio of 
factorials, reference [102] gives the conditions of these factorials for the hypergeometric 
series to be algebraic (all the coefficients are thus integers). A simple example is, for 
instance the algebraic function: 

:>8 



3F2 



n=0 



Along this line it is worth recalling Delaygue's Thesis [93] (see also 

Bobcr [104]) which gives some results^ for series expansions^ such that 

their coefficients are ratio of factorials, namely 2^1 ([1/3, 2/3], [1]; 27 x) , and 
4 F 3 ([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; 256 x), giving respectively the series 

^ (3n)! „ ^ ((2n)!) 4 „ 

2^ 7^)3 ' x > 2^ -Trip- ■ x > and: 

n=0 v 1 n=0 y ' 



These ratio of factorials are all integer numbers. 



Appendix C. Yukawa couplings 

Appendix C.l. Yukawa couplings as ratio of determinants 



Consider an order-four MUM linear differential operator. Let us introduce the 
detcrminantal variables W m — det(M m ) which are the determinants! of the following 
m x m matrices M m , m = 1, . . . , 4, with entries expressed in terms of derivatives of 
the four solutions yo(x), y\{x), y2{x) and 2/3(2;) of the MUM linear differential operator 



(see Section (4.1) for the definitions). One takes W\{x) 



Mo 



where: 



2/o 
2/o 



2/1 

y'i 



Ah 



dx ' 



2/0 
2/o 

y'o 



2/1 
y[ 
y'i 



2/2 
2/2 
2/2' 

da; 2 ' 



yo(x) 
yo 



and: 
2/1 
2/i 
2/i' 



2/2 
2/2 
//2 



Since q is equal to q = exp(y 1 /yg), and its derivative verifies 



2/i- 



^ dg 



1U 2 



d 

c/.r 



w 2 



d_ 

dx ' 



we have that 



q ■ 



dq 



dx 2 



W| dx 



d_ 

dx 



dW 2 
dx 



d 

dx 



2/3 
2/3 

2/3' 



2/b 2/i 2/2 2/3 
d 3 



(C.l) 



(C.2) 



(C.3) 



(I Necessary and sufficient conditions for the integrality of the mirror maps series. 

1 These series are not algebraic functions. 

X For an order- four operator the Wronskian is W4. 



K( n \ ^(g) 3 • w l' 2 ( 1 Mq) \ 3 Wl /2 ( dx(g) y 



(C.6) 
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We deduce, after some simple algebra, an alternative definition for the Yukawa 
coupling: 

K(g) = = ^ = (C4) 

W V dq) \y Q J W$ W$ V ; 

to be compared with the other previous alternative expression previously given (86) 
for the Yukawa coupling 

) 3 

Vl ^x(q) dq J yl V dq ) 

In fact from (C.2) we deduce 

^(g) \ 3 w[ y$_ 

dq ) W$ Wi' 

and, so, (C.5) is compatible with (C.4) if the following identity is verified: 

Wl = Wi-yl = Wi- Wl. (C.7) 

This identity is in fact specific of order-four operators conjugated to their adjoints (see 
below (C.17)). Therefore we prefer to use definition (C.4) for the Yukawa coupling, 
instead of the more restricted definition (C.5). 

Let us assume that the pullback p(x) has a series expansion of the form 

p(x) = A- x r ■ A(x), (C.8) 

where the exponent r is an integer, where A is a constant, and where A(x) is a 
function analytic at x = with the series expansion: 

A{x) = 1 + o?i • x + a 2 ■ x 2 + ■ ■ ■ 

The detcrminantal variables W m transform very nicely under pullbacks p{x) of the 
form (C.8): 

(W 1 (x), W 2 (x), W 3 (x), Wi{x)) — > (C.9) 

(Vi(p(z)), - ■ W 2 (p(x)), ^ • W 3 ( P (x)), ^ ' W 4 (p(x))), P ' = 

V r r° r° / ax 

One can show that the nome (36) of an order- N operator transforms under a 

pullback p{x): 

q{x) — > Q(x) with: X-Q(x) r = q(p(x)). (CIO) 

From the covariance property (C.9), and from the previous transformation 
q — > A • q r for the nome, one easily gets the transformation of the Yukawa coupling 
seen as a function of the nome K(q) — > K(X ■ q r ): 

Wx jxf ■ W 3 (x) 
W 2 {xf 

w^p^y -w 3 {p{x)) 

W 2 {p{x)f = K ( q W x W = K(X-Q(x)). 

For A = 1 and r = 1 (i.e. when the pullback is a deformation of the 
identity transformation), one recovers the known invariance of the Yukawa coupling 
by pullbacks (see Proposition 3 in [105]). 

One finds another pullback invariant ratio, namely: 



K{q{x)) = (C.ii) 



Wi ■ W$ 

K = 1 3 _ , (C.12) 
W 4 ■ W? ' 
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which is, in fact, nothing but the Yukawa coupling for the adjoint of the original 
operator. 

Another invariance property is worth noting. Let us consider two linear 
differential operators f2i and Q 2 of order N that are equivalent, in the sense of 
the equivalence of linear differential operators. This means that there exists linear 
differential operators intertwincrs I\, I 2 , Ji, J 2 , of order at most N — 1 such that 

Ql ■ h = h- fl 2 , and: Ji • O x = fi 2 • J 2 . (C.13) 

Let us assume that one of these intertwincrs is a linear differential operator of order 
zero (a function), then the previous homomorphism between operators amounts to 
saying that the two operators are conjugated by a function: 

2 = p(x) ■ Hi • p(x)~\ (C.14) 

which correspond to changing the four solutions as follows: yi — ¥ p(x) ■ j/j. In such 
a case the previous determinant variables transform, again, very nicely under the 
"gauge" function p{x): 

(Wi, W 2 , W 3 , W 4 ) ->■ (p(x) ■ Wi, p(x) 2 ■ W 2 , P {xf ■ W 3 , P {xf ■ Wa). (C.15) 

It is straightforward to see that the Yukawa coupling and the "dual Yukawa", 
K* , are invariant by such a transformation^. Two conjugated operators (C.14) 
automatically have the same Yukawa coupling. 

Do note that the Yukawa couplings for two operators, which are non trivially 
homomorphic to each other (intertwincrs of order one, two, ...), are actually different. 
The (pullback invariant) Yukawa coupling is not preserved by operator equivalence (see 
subsection (8.2)). 

Remark: The definition of these determinantal variables Wi heavily relies on 
the MUM structure of the operator between the four solutions the definition of WiS, 
in particular the log-ordering of the solutions. It is worth noting that if one permutes 
the four solution j/j, one would get 24 other sets of (Wi, W 2 , W3, W4) which are 
actually also nicely covariant by pullbacks, thus yielding a finite set of other "Yukawa 
couplings" or adjoint Yukawa coupling K* also invariant by pullbacks. 

In fact these "Yukawa couplings" (C.4), and other adjoint Yukawa K* (C.12), 
can even be defined when the linear differential operator is not MUM, and they are 
still invariant by pullbacks. 



Appendix C.2. Fullback-invariants for higher order ODEs 

These simple calculations can straightforwardly be generalised to higher order linear 
differential equations. We give here the invariants for higher order linear differential 
operators. 

Let us give, for the n-th order linear differential operator the list of the K n 
invariants by pullback transformations: 

wi 1 ■ W 3 W* ■ Wj wl 5 ■ w 5 W 2A ■ w 6 

3 ~ wi ' 4 w$ ' 5 ~ W}° ' 6 " wi 5 ' " ' 

K n = W \l' b Wn > with: a « = n-(n-2), b n = n < n ~ l \ (C .16) 
W 2 n 1 

*i K and K* (and their combinations) are the only monomials W" 1 W£ 2 W^ 3 W^ 4 to be invariant 
by (C.9) and (C.15). 
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A n-th order linear differential operator has K n as an invariant by pullback 
transformation, as well as all the K m with m < n. K s is the Yukawa coupling, 
and one remarks, for the order-four operators, that the other pullback invariant K* 
(see (C.12)), which is actually also the Yukawa coupling of the adjoint operator, is 
nothing but A'f/A 4 . 

For order-four operators conjugated to their adjoint (see (C.14)) (i.e. operators 
homomorphic to their adjoint, the intertwiner being an order zero differential operator, 
a function), one has the equality 

Ki = Kl, i.e. K 3 = K* , or W$ = W 2 ■ W 4 , (C.17) 

to be compared with the equality in Almkvist et al. (see Proposition 2 in [70]) 

2/02/3 - ysy'o = yiy' 2 - 2/22/1, (C.18) 

which is satisfied when the Calabi-Yau condition that the exterior square is of order 
five is satisfied. 

If a linear differential operator verifies condition (C.17), its conjugate by a 
function, p(x) ■ ^4 • p(x)^ 1 , also verifies condition (C.17) (their Yukawa couplings are 
equal) . 

The condition (C.17) is not satisfied for linear differential operators homomorphic 
to their adjoint with non-trivial intertwiner (of order greater than zero). For instance 
the order-four operator (87) does not satisfy condition (C.17). 

Appendix D. More Hadamard products: a Batyrev and van Straten 
Calabi-Yau ODE 

An order-four operator has been found by Batyrev and van Straten [72] associated 
with a Calabi-Yau three-fold on P2 x P 2 

Bi = 6 4 -~3x-(76 2 + 79 + 2) -{30 + l)-(3 + 2) (D.l) 

- 72 x 2 ■ (3 6* + 5) • (3 6* + 4) • (3 6* + 2) ■ (3 6* + 1). 

This operator is conjugated to its adjoint: B\ ■ x = x ■ adjoint(Bi). 

Operator (D.l) is a Calabi-Yau operator [72]: it is MUM, and it is such that its 
exterior square is of order five. It has a solution analytical at x = which is actually 
the Hadamard product of the previous selected hypergeometric 2F1: 

^([i, §],[!]; 27.) * • 2 F,([I, |], [1]; -J^)). (D.2) 

The coefficients of the series expansion of (D.2) are integers: they can actually be 
written as a sum of product of binomials: 




(D.3) 



The series expansions of the nome and the Yukawa coupling of this Calabi-Yau 
operator are series with integer coefficients, reading respectively: 

q = x +48 x 2 + 4626a; 3 + 549304 x 4 + 74589735s 5 + 11014152048 a: 6 + ••• 
K(q) = 1 + 21 q + 3861 q 2 + 429159 q 3 + 57224661 q 4 + 7337893896 q 5 + • • • 
K{x) = 1 + 21 x +4869 x 2 + 896961 x 3 + 175176657 x 4 + 34770008997 x 5 + ■■■ 
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